Electrostatic potential profiles of molecular conductors 



G. C. LiangB A. W. GhoshQ M. Paulsson, and S. Datta 
School of Electrical and Computer Engineering, Purdue University, W. Lafayette, IN 4^907 

(Dated: February 2, 2008) 

The electrostatic potential across a short ballistic molecular conductor depends sensitively on the 
geometry of its environment, and can affect its conduction significantly by influencing its energy 
levels and wave functions. We illustrate some of the issues involved by evaluating the potential 
profiles for a conducting gold wire and an aromatic phenyl dithiol molecule in various geometries. 
The potential profile is obtained by solving Poisson's equation with boundary conditions set by the 
contact electrochemical potentials and coupling the result self-consistently with a nonequilibrium 
Green's function (NEGF) formulation of transport. The overall shape of the potential profile (ramp 
vs. fiat) depends on the feasibility of transverse screening of electric fields. Accordingly, the screening 
is better for a thick wire, a multiwalled nanotube or a close-packed self-assembled monolayer (SAM), 
in comparison to a thin wire, a singlewalled nanotube or an isolated molecular conductor. The 
electrostatic potential further governs the alignment or misalignment of intramolecular levels, which 
can strongly infiuence the molecular I-V characteristic. An external gate voltage can modify the 
overall potential profile, changing the current-voltage (TV) characteristic from a resonant conducting 
to a saturating one. The degree of saturation and gate modulation depends on the availability of 
metal-induced-gap states (MIGS) and on the electrostatic gate control parameter set by the ratio 
of the gate oxide thickness to the channel length. 

PACS numbers: PACS numbers: 85.65.+h, 73.23.-b,31.15.Ar 
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Recently there has been considerable progress in the 
experimental analysis Q, IE Q and theoretical model- 
ing 0, IE S H of molecular electronic devices. Con- 
duction through a molecule depends on its intrinsic chem- 
istry, as well as external influences such as the electrode 
geometry, charging and bonding at the electrodes. In 
particular, the electrostatic potential profile in molec- 
ular devices is a quantity of great interest, because it 
contributes to the self-consistent field in the molecular 
Hamiltonian, thereby influencing the electronic proper- 
ties and regulating the flow of current. The spatial vari- 
ation of the electrostatic potential carries nontrivial in- 
formation about molecular screening, the presence of im- 
purities and Schottky barriers, as well as features related 
to the alignment of energy levels within the molecule 
[lolllll |. Furthermore, the long-ranged nature of electro- 
static forces allows the molecular levels and wavefunc- 
tions to be tuned remotely, with a gate electrode for ex- 
ample. 

It is often stated that the voltage drop in a ballistic 
conductor needs no discussion, for all the drop must be 
localized at the interfaces with the contacts. While this 
is true in some sense for the electrochemical potential pro- 
file, fi{r), it is certainly not true for the electrostatic po- 
tential, (j^ir); in fact, the two profiles can differ widely 
in low-dimensional materials with long screening lengths. 
We will try to bring out this distinction later using a sim- 
ple classical diffusion model, since fj,{r) is a tricky concept 
to define under more general conditions of transport. The 
electrostatic potential ^(r^, by contrast, is a clearly de- 
fined concept even for quantum transport far from equi- 
librium. This paper is essentially about (/)(r). Under 
certain conditions the electrostatic potential profile can 



have significant effects on the current-voltage character- 
istics of a conductor ^El- In ^ standard two-terminal 
molecular conductance measurement for example, 

the precise division of the applied bias between the source 
and drain contacts can cause the I-V to change from sym- 
metric to asymmetric, with the conductance gap deter- 
mined either by the molecular levels exclusively or by the 
contact Fermi energy in addition |l3l |. Such a sensitive 
dependence of current conduction on the potential pro- 
file is not typical in mesoscopic physics, but quite routine 
in molecular electronics. We will try to provide simple 
insights to describe such conditions. 

Direct measurement of the electrostatic potential pro- 
file of a molecular conductor is challenging, given its 
small size. Attempts at direct AFM-based or potentio- 
metric measurement of the profile have been limited to 
long (~ 0.3 — 3 um) carbon nanotubes or organic 
molecular solids where the profile has been deter- 
mined mainly to be flat inside the molecular system, with 
the voltage drop largely occurring at the contacts. Simi- 
lar profiles have been theoretically postulated or invoked 
in various semi-empirical treatments of molecular con- 
duction IE ll^ . In sharp contrast, a ramp- like poten- 
tial across a molecular wire has been calculated by semi- 
empirical ^E| J as well as a number of first-principles den- 
sity functional theory (DFT)-based simulations |E S ■ 
It is further reasonable to expect that potential variations 
on an atomic scale would be influenced by atomistic fea- 
tures on the contacts and surrounding molecules. The 
nature of the potential proflle thus needs to be sorted 
out. 

In this paper we perform a fully quantum kinetic, 
atomistic treatment of the electrostatic potential proflle 
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FIG. 1: Schematic of two-terminal geometry, with metal con- 
tacts of infinite cross section coupled to an active device. 
The device will be replaced by a gold wire of various shapes 
throughout the rest of the paper, and by a phenyl dithiol 
(PDT) molecule towards the end. Part of the contact mate- 
rial may be included in the device in addition. 

across a prototype molecular conductor, and examine its 
influence on experimentally measurable current-voltage 
(I-V) characteristics of the conductor. The prototype 
conductor we investigate is either a gold atomic chain 
or a phenyl dithiol molecule, attached to a cluster of 
metal atoms on both ends and sandwiched between metal 
contacts of infinite cross section (Fig. The "trivial" 
Laplace part (j)L of the potential profile depends on the 
geometry and dielectric constants of the molecule and 
the electrodes. This part is expected to dominate the 
electronic properties for insulators and molecular con- 
ductors away from resonance, and can be modulated ex- 
ternally with a gate, leading potentially to transistor 
action. Furthermore, the Laplace part determines the 
alignment of energy levels localized on different parts of 
the molecule, which could drive interesting quantum ef- 
fects such as current rectification or negative differential 
resistance (NDR). In contrast, the quantum contribu- 
tions to the potential through the Poisson term (t>p{p) 
are affected by the atomicity of the device. The Poisson 
part becomes important in conducting systems such as 
metallic wires, metallic nanotubes and molecules driven 
into resonance with a source-drain bias or turned on with 
a gate. 4>p{p) contains nontrivial physics, including self- 
consistent charging, screening and Friedel oscillations. 

Section I provides a semiclassical description of con- 
duction that outlines the difference between the electro- 
chemical and the electrostatic potentials. Section II de- 
scribes the nonequilibrium transport formalism that we 
employ, and the model geometries and Hamiltonians that 
we use. Section III calculates the voltage drop within the 
contacts, illustrating that only a fraction of the applied 
bias drops across a charge neutral region including the 
molecule and the contacts. Section IV discusses the de- 
pendence of the potential profile on the thickness of the 
wire, specifically addressing the feasibility of adequate 
amount of transverse screening. Similar issues related 
to ineffective screening in 1-D have been raised in the 
context of carbon nanotubes 0|. Section V addresses 
the importance of screening by the environment, such as 
in a self-assembled monolayer (SAM) of molecular wires. 



In section VI we look at potential variations within the 
molecule, which although small could nevertheless signif- 
icantly influence observable electronic properties, such as 
the current-voltage (I-V) characteristics. Finally in Sec- 
tion VII, we talk about the influence of remote contacts 
such as a gate electrode on the potential profile, which 
in turn determines three-terminal transistor I- Vs. We 
summarize our results in section VIII. 

I. SEMICLASSICAL DESCRIPTION OF 
POTENTIAL PROFILE 

In this section, we present a simplified semiclassical, 
continuum description of current conduction, invoking 
macroscopic parameters such as the conductivity tensor 
and the dielectric constant. The aim of this section is 
to provide an elementary description of the distinction 
between electrochemical and electrostatic potentials. For 
molecular systems one cannot use such a macroscopic 
description, and a proper quantum kinetic description 
needs to be invoked (section II). 

A semiclassical description involves defining a conduc- 
tivity tensor cr(r). The electrochemical potential satisfies 
the equation of continuity, while the electrostatic poten- 
tial satisfies Poisson's equation: 

V • (aVfi^ = (la) 

V • (eV^) = [n{r) - No] (lb) 

where e is the dielectric constant. No is the dopant den- 
sity, and n(r) is the electron density: 

n(r) = j dED{E, r)h{E + 4>{r) - ti{r)) 

UE) = [l + cxp(i?/fcBT)]-\ (2) 

D{E,r) representing the local electron density of states 
(LDOS). A conductivity mismatch at the device-contact 
interfaces, such as that generated by a variation in the 
cross-sectional geometry or the doping profile, allows us 
to hold the contact electrochemical potentials at fixed 
voltages under bias, dropping u almost entirely at the 
interfaces for a ballistic device [l3| • This variation (5/x in- 
fluences the local charge distribution 5n{r) through the 
chemical potential /x— </> (Eq.[5J), so that the electrochemi- 
cal potential profile ^(r) in effect acts as the driving force 
for the electrostatic potential profile 5(j) in Eq. Ilbl In 
linear response, 5n{f) « Dq {5^, — S(f>), where Dq is the 
density of states (DOS) at the Fermi energy. Poisson's 
equation ljlb|l then reads: 

(v^-^).0=-^^, (3) 

indicating that S(j) is given by a convolution of Sj-i and 
a screening function that varies on a Dcbye (Thomas- 
Fermi) lengthscale given by Ad ~ -^/e/e^Dg Um- (^i^ 
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analogous expression can be invoked to describe screen- 
ing by surface states |2Cll|'). The electrostatic potential thus 
has a slower spatial variation than the electrochemical po- 
tential profile \2l\ . In metallic conductors with a high 
DOS and correspondingly small Debye length the two 
profiles track each other in order to avoid large charge 
buildups. In contrast in a semiconductor or insulator 
having small DOS inside the bandgap the Debye length 
can be quite large, so that the two potentials can vary 
on widely different length scales A molecular con- 

ductor is intermediate between the two limits, acting as 
an insulator when the contact electrochemical potentials 
are off-resonant with the levels, and as a conductor on- 
resonance. One can thus have a nontrivial electrostatic 
potential profile even in a ballistic molecular device where 
the electrochemical potential does not vary spatially at 
all, except at the contact-molecular interfaces. 

Having illustrated the basic distinction between the 
two kinds of voltage-drop, we now move onto a rigorous 
quantum kinetic and atomistic description of it. 



matrix p{f, r*) of which the electron densities na{r) form 
just the diagonal parts. A formal way of handling these 
issues is by using the nonequilibrium Green's function 
(NEGF) prescription 

In the NEGF formalism, one deals with retarded, ad- 
vanced, lesser and greater single-particle Green's func- 
tions [lal • The retarded Green's function matrix is 

GiE) ^[ES-H- Uscfip) - ^siE) - ^oiE)]-' . (5) 

where H is the single-particle molecular Hamiltonian in 
an appropriate basis set, S is the overlap matrix in that 
basis-set, the self-consistent potential Usct is dominated 
by the electrostatic potential </> (and in principle, includes 
exchange-correlation effects too), and the self-energy ma- 
trices T,s,D represent the influence of scattering by the 
source (S) and drain (D) contacts. The contact self- 
energies can be calculated for electrodes of given geom- 
etry and surface bonding 0, yielding level broadenings 



II: QUANTUM KINETIC FORMALISM 



S.D 



^S.D 



(6) 



Basic equations. Eq. Hal needs to be modified to solve 
for the full quantum transport under bias. At equilib- 
rium, one could still use the concept of an electrochemical 
potential, defined by the contact Fermi energy Ep. The 
electronic local density of states (LDOS) are obtained by 
solving Schrodinger's equation for the molecular Hamil- 
tonian, supplemented by a self-consistent potential C/scf, 
and filling the eigenstates according to equilibrium sta- 
tistical mechanics: 

[H + Uscf] ^-olf) = S„*„(f) 

"cq(r) = ^n„(f) = l*«(r)P/o(i?„ - Ep). (4) 

a a 

In general, the self-consistent potential C/scf consists of 
a Hartree-term obtained from the electrostatic potential 
(j> that is the solution to Poisson's eauation llbl with n re- 
placed by ricq, and additional exchange-correlation con- 
tributions that need to be incorporated through an ap- 
propriate ab-initio technique. The expression for n can 
be recast in the form Eq.Elby recognizing that the LDOS 
DiE,r)=j:j^.m'SiE-E^). 

Under bias, the contact electrochemical potentials sep- 
arate, and the system is driven out of equilibrium in its 
bid to equilibrate with both contacts, causing a current 
flow. Under these nonequilibrium conditions an unam- 
biguous common electrochemical potential is hard to de- 
fine. One could, however, describe transport in terms 
of groups of electrons that are separately in equilibrium 
with the two contacts Furthermore, the contacts 

could fill the energy levels in a correlated way, so that 
one needs to deal with the full nonequilibrium density 



The nonequilibrium density matrix p determining 0(p) is 
calculated self-consistently within the NEGF formalism 
using the lesser Green's function as 0, 0| : 



j dE [-iG<[E) /2tt] 



(7) 



with the electron density representing its diagonal com- 
ponent (see Eg. If 2|l . The function —iG'^{E) describes 
how the molecular states are filled in a correlated way by 
the two contacts, and is itself determined by the retarded 
Green's function G (Eq. |SJl, the broadenings Vs.d-, and 
the contact electrochemical potentials ps.D- 



fsAE) = 1 + cxp 



G [fsTs + fnTD] Gt 
E - ^IS.D 



(8) 



Since we explicitly include any asymmetry in electrostatic 
coupling with the contacts through (/'(p) and in quantum 
coupling through Ts,d, we dispose ^s,d symmetrically 
about the equilibrium contact Fermi energy Ep under 
a source-drain bias V: p^s.D = Ep =p eV/2. The con- 
tact electrochemical potentials help fix the boundary con- 
dition (Dirichlet or Neumann) for Poisson's equation in 
the contacts, so that the electrostatic potential (t>{p) within 
the device can now develop on its own self-consistently 
under the applied bias |2,ll |. 

The NEGF density matrix equations CEas. iyiHl) provide 
the nonequilibrium generalization of Eq. 0] The corre- 
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spondence is evident by rewriting Eq. ^ as 



[Peq] 



J dEfo{E - Ef)D{E) 



dE 
2^ 



(9) 



where [p] and D are the charge density matrix and the 
density of states matrix whose diagonal entries give the 
LDOS and the trace gives the total DOS. The self- 
consistently converged Green's function G is then used 
to obtain the current. For coherent transport the NEGF 
current expression formally resembles that in Landauer 
theory ^ . 

2e 

I = dET{E)[fs{E)- fD{E)] 

T = Trace (rsGTcGt) . (10) 

Given an appropriate Hamiltonian iJ, self-consistent po- 
tential (f) and self-energies S, eans. l5l6l7l8ll0l vicld a com- 
plete set of equations allowing us to obtain the nonequi- 
librium charge and current densities under bias. Al- 
though the equations describe coherent transport, NEGF 
allows us to handle incoherent processes through addi- 
tional self-energy matrices determined by the correspond- 
ing scattering potentials. 

Incoherent processes. Incoherent processes such as 
hopping or phonon scattering can be included rigorously 
through a self-energy in the NEGF prescription 01 • A 
simple, phenomenological way of including this would 
be to model each scattering center 0jS 0, floating voltage 
probe with which the electron equilibrates locally, moti- 
vated by Biittiker [isl [isj . Such a probe docs not draw 
any net current, but randomizes the phase of the incom- 
ing electrons by reinjecting them into the device. The 
probe can be parametrized by two quantities, its elec- 
trochemical potential and the self-energy Ep. The 
Green's function expression in Eq. [S] includes the addi- 
tional self-energy contribution Ep, while the transmis- 
sions Tij between any of the three contacts is given by 
T,j = Trace (r^GF^Gt), i,j^S,D,p. 

Besides modifying the transmission, incoherent scat- 
tering also affects the local electron density through /ip, 
which is fixed by requiring that the net current drawn 
by the voltage probe is zero. One could model phase- 
breaking without energy relaxation^ where the probe 
Fermi functions fp{E) are adjusted at each energy such 
that the current Ip{E) drawn by the probe at every in- 
dividual energy channel is zero. The expression for the 
net transmission in this case becomes relatively straight- 
forward: 



T{E)^TsDiE) 



Tsp{E)TpD{E) 
Tsp{E)+TpD{E) 



(11) 



Alternately one could include energy-relaxation, adjust- 
ing the overall probe electrochemical potential such 
that the net probe current Ip = J dEIp{E) obtained 
by integrating the current contributions over various en- 
ergy channels adds up to zero. One will then need to 
modify Eq. ^| to include the contributions from the 
probe to the net source-drain current. The probe self- 
energy ^p{E) can be calculated using various models for 
scattering relaxation. For phonons for example, we use 
the self-consistent Born approximation whereby the 
phonon self-energy depends recursively on the electron 
Green's function. In the limit of low-energy phonons, 
this amounts to setting Ep = AqG, Aq being the cou- 
pling between the molecule and the phonon bath. In 
gen eral, however, the expressions are more complicated 

Model Hamiltonians for device and contacts. The 
above formalism is general, and requires an appropriate 
Hamiltonian H describing the intrinsic molecular chem- 
istry, an adequate treatment of the geometry and bond- 
ing at the contact surfaces described by self-energy ma- 
trices J^s,D, as well as a suitable self-consistent poten- 
tial describing the electrostatics of the device. In the 
past, we calculated these matrices by modifying a stan- 
dard quantum chemical software (GAUSSIAN98), self- 
consistently coupled with a NEGF-based transport for- 
malism 0. Such a modification allowed us to get I- 
V characteristics and potential profiles using a density 
fimctional (DFT) description of both the molecule and 
the contacts. While DFT codes allow 'first principles' 
treatment of electron-electron interactions and quantum 
correlations with no adjustable parameters, they are 
tedious and time-consuming. Furthermore incorporat- 
ing the Poisson boundary conditions typically requires 
a large chunk of the contacts included along with the 
molecule within our device |26l | . Alternatively, one needs 
to perform repeated transformations between an orthog- 
onal real space basis where Poisson boundary conditions 
are readily incorporated, and a non-orthogonal orbital 
basis suitable for describing the atomic bondings and the 
molecular chemistry. 

In order to bypass the complexities associated with 
implementing boundary conditions within DFT, we will 
concentrate here on a much simpler semi-empirical de- 
scription of the device and contacts, described by an or- 
thogonal basis set {^^} of one orbital per atom. Such 
a basis could, for instance, describe the essential physics 
of molecular transport through conjugated Pz (tt) elec- 
tron systems, and s electrons for metallic gold wires. In 
this paper, the scalar on-site and hopping parameters are 
chosen to match experimental values of the Fermi energy 
and density of states for gold, as well as the energy eigen- 
values and orbital shapes for phenyl dithiol (PDT) calcu- 
lated u sing DFT within the local density approximation 
(LDA) |2a (parameters listed in TablcHl. The source and 
drain contacts arc modeled as Au(lll) surfaces, with the 
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On-site 
(eV) 


Hopping 
(eV) 


Hubbard 
(eV) 


Bond-length 
(A) 


Au-Au 


-4.3 


8.75 


11.13 


2.885 


C-C 


-4.42 


2.5 


11.13 


1.4625 


S 


-6.32 




9.94 




S-C 




1.5 




1.8 


S-Au 




1.6 




2.885 



TABLE I: Semi-empirical parameters used to simulate molec- 
ular wires. The S-Au coupling is reduced from 1.6 to 0.8 
in the last section to simulate gating of a weakly contacted 
molecular wire. 



wire bonded symmetrically to a surface triangle of gold 
atoms. The contact self energies Ss,d are obtained rig- 
orously us ing a real space recursive formalism described 
elsewhere lEllIllil. 

Self- consistent potential. The self-consistent potential 
4>{p) in Eq. [S] can be obtained for a given charge density 
matrix p in a variety of ways. We will ignore the evolu- 
tion of the exchange-correlation potential under bias (the 
equilibrium result is incorporated in H) , and concentrate 
only on the Hartree part |29j , which depends only on the 
diagonal elements n of the density matrix. A computa- 
tionally tedious but accurate way to obtain the Hartree 
contribution is the direct solution of 3-D Poisson's equa- 
tion numerically on a real space grid, with appropriate 
Dirichlct/Neumann boundary conditions at the bound- 
aries (section IV). A faster way, that involves grid points 
only on the charges and the contacts is the method of 
moments (MOM) [s^l • To understand this, one starts by 
writing the solution to Poisson's equation in matrix form 
in a suitable basis {^}, and then partitioning the system 
into the device ('d') and the contact ('c') sectors: 



4>d 



Udd Udc 

Ucd Ucc 



nd 
ric 



n{r) = p(f, T) - %i^P^.'^u{r^^ (12) 

U representing the Hartree term [sJl- Eliminating the 
contact charge density ric, one gets the device electro- 
static potential that is used in Eq.|31by solving the matrix 
equation: 

Hp) = ^d = UacUcc^^c+{Udd-UdcU-c'U,a)n^ 

= cPl + Mp), (13) 

where the Laplace part is the solution to the homoge- 
neous Poisson's equation corresponding to zero charge 
in the device (rid = 0) and specified contact potentials, 
while the inhomogeneous Poisson part corresponds to the 
solution for given contact and device charges and with 
zero potential on all the contacts {(j)c — 0). The advan- 
tage of MOM is that we need grid points only on the 
contacts with specified potentials 4'c and on the device 
with specified free charges n^. The formalism can easily 



be extended to include polarization charges in dielectric 
materials. 

In this article, the Coulomb matrices U are obtained 
using the Pariser-Parr-Pople (PPP) model within the 
Matago-Nishimoto approximation |32l |: 



1/7. 



(14) 



where 7,,, = Q.5{Uii + Ujj)/e'^. The above term interpo- 
lates between asymptotic 1/r Coulomb repulsion at large 
distances and on-site Hubbard repulsion Ua at short dis- 
tances (table nj (Alternate ways of handling the Hub- 
bard term exist; see for e.g. Mc. Lennan et al., [2]| . 
Appendix). The effect of image charges on the con- 
tact surface atoms is already included in Eg. II 31 through 
the second term in brackets in the expression for 4>p{p). 
Additional, possibly negligible, image contributions from 
deeper within the contacts are added in by hand by gen- 
erating a series of contributions similar to Eq. 1141 due to 
electrostatically infinite planar contacts • 

The boundary conditions on Poisson's equation are 
set by the electrochemical potentials p-s.D- One can ei- 
ther use Neumann boundary conditions, requiring overall 
charge neutrality deep inside the contacts (section IV). 
Alternatively, one could enforce Dirichlet boundary con- 
ditions, setting the electrostatic potentials equal to the 
electrochemical potentials deep inside the contact (rest 
of the paper). 

Having laid out our general framework and the de- 
scription of our model Hamiltonians, we will first try to 
determine how much of the applied bias across a metal- 
molecule-metal heterostructure drops within the contacts 
themselves. 



III. HOW MUCH VOLTAGE DROPS ACROSS 
THE MOLECULE? 



The transport formalism described in section II applies 
equally to the device and the contacts. To distinguish 
between the two, one needs a dilution of modes in going 
from the contact to the device, such as a drastic change 
in doping profile (Figs. |21i,b) or a geometrical narrowing 
of the cross-section (Fig. EJ;). A slight deviation from 
equilibrium inside the contact, integrated over its multi- 
ple modes, is then sufficient to drive current through the 
device, while keeping the contact electrochemical poten- 
tials essentially constant. In terms of Eq. ^| one needs 
a large conductivity inside the contacts relative to the 
molecule in order to hold the electrochemical potentials 
constant inside each contact and drop the applied volt- 
age bias across the molecule itself. For a ballistic device, 
however, part of the applied voltage drops within the 
contacts, as we now describe. 

To simulate the potential profile under bias across the 
device as well as the contacts, we will consider an ex- 



6 




FIG. 2: Electrostatic potential profile in various structures under zero-field boundary conditions (ZFB). Part of the contacts 
(blue) has been incorporated into the device region within dashed lines, (a) Charge neutrality restrictions due to ZFB in a 
silicon (n^nn^) structure causes the source electrostatic and electrochemical potentials under a 0.2 V bias to separate, thus 
dropping some of the applied voltage in the contacts. The voltage "loss" can be reduced by (b) increasing scattering or by (c) a 
geometric dilution (see text), which reduces this mismatch. The structures in (a) and (b) have the same geometry and doping 
properties, shown in the figure. Structure (c) is simulated by using gold parameters (Table 1). The channel consists of a four 
atom gold chain, while each contact cluster within the device region consists of a 10{Lx) x ll(Ly) simple cubic gold lattice. 



tended device that actually incorporates a large enough 
part of the contacts to keep it essentially charge neu- 
tral. Such an inclusion would automatically take care of 
image charges for metal contacts or depletion/inversion 
charges in semiconducting contacts, for example. How- 
ever, charge neutrality ensures that only a fraction of the 
applied bias appears across the ends of the device, es- 
pecially if it is ballistic, the rest dropping deeper inside 
the contacts. It is hard to impose Dirichlet boundary 
conditions on the electrostatic potential at the two ends 
without knowing this "lost" fraction a-priori; a better 
way seems to be to just fix the electrochemical potentials 
fJ-s,D cit the two ends, impose Neumann-type zero-field 
boundary conditions (ZFB) on the electrostatic potential 
consistent with Gauss ' law and charge neutrality, and let 
the potential profile develop on its own self- consistently. 

To illustrate the voltage loss within the contacts de- 
scribed above, we simulate a n-doped device (Fig. |21[a)) 
consisting of a cluster of silicon atoms of infinite cross- 
section attached to a n-|--|- silicon contact of identical 
structure and cross-section but with a heavier doping 
density. Part of the contact has been incorporated into 
the extended device simulation region (dashed lines) to 
account for surface charges. The equilibrium {V = 0) 



potential profile shows a barrier on the device part that 
reflects the incoming electron waves from the contacts. 
We determine the potential profile under bias by explic- 
itly solving 1-D Poisson's equation self-consistently on a 
real space grid using the Newton-Raphson method [s^l 
with ZFB at the ends. Note also that in the transport 
cauations l8l9ll0l onc needs to replace each Fermi function 
f{E) by the corresponding 2-D Fermi function F2Ti{E) 
obtained by summing the transverse modes over the two 
irrelevant dimensions [33,[33|. As is evident from Fig.|3i, 
applying a 0.2 V source-drain bias lowers the drain elec- 
trochemical potential, which in turn lowers the electro- 
static potential barrier and leads to current flow. The 
net self-consistent potential, (j){V = 0.2) — 0(F = 0) 
across the wire is shown in solid line, and exhibits signif- 
icant screening due to charge rearrangement within the 
n-dopcd region. Significantly, lowering the drain elec- 
trochemical potential also lowers the source electrostatic 
potential, implying that only a fraction of the applied 
0.2 V bias actually appears across the simulated region. 
To understand this, we recall that in a ballistic device 
the current carrying electrons can be grouped into -l-fc 
and —k states that originate in the drain and source con- 
tacts respectively 0|. A bias changes the occupancies 
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of these states, so that a current flows due to their im- 
balance. Inside the source contact too, a bias produces 
a similar imbalance between +k and —k states, with the 
result that part of the electron states inside the source 
end up being depopulated by the drain electrode. The 
depopulation forces part of the applied bias to drop in- 
side the source contact in order to keep it charge neutral 
|37l |. Consequently, only a fraction of the applied bias 
occurs across the two ends. 

The depopulation of the source contact is particularly 
noticeable for ballistic transport, where the drain end can 
effectively empty out the source end and lower its electro- 
static potential significantly. This observation seems to 
be in keeping with Landauer's idea that only a fraction 
R {R: device reflection coefficient) of the applied bias 
actually appears across the active device, which in con- 
junction with a current proportional to the transmission 
T gives the celebrated Landauer conductance for a single- 
moded structure minus the interfaces, G = (2e^//i)(T/i?) 
[ssl l . It seems evident therefore that the way to eliminate 
the voltage loss for a ballistic device is to increase the 
average scattering (i? ^ 1), which would in effect isolate 
the source and drain and prevent the latter from depop- 
ulating the former. This could be achieved by increasing 
the doping density in the n++ region of Fig. I3a), by 
increasing scattering throughout the device (Fig. I^^b)), 
or by enforcing a geometric dilution with a flared-out 
contact geometry (Fig. Hfc)). In (b) we use Biittiker 
probes to simulate scattering, while in (c) we solve the 2- 
D Poisson equation to take care of the varying transverse 
cross-sectional geometry, using the relevant -Fid function 
obtained by summing the Fermi-Dirac function over the 
single irrelevant dimension. As we see in both cases, a 
larger fraction of the applied bias appears across the ends 
with a smaller voltage loss in the contacts, as evidenced 
by the closer agreement between the electrostatic poten- 
tial values at y = and V = 0.2 V deep within the 
source [s^. Scattering in the contacts can thus have a 
significant influence on the device conductivity, tending 
to make it robust with respect to spatial variations in the 
intcrfacial geometry [4fil |. 

In the simulations we present hereafter, we will ex- 
plicitly incorporate a built-up, atomistic part of the con- 
tacts along with the molecule within our device, enforc- 
ing thereby a geometric dilution. Accordingly, we will 
assume that most of the applied bias appears across 
the ends, with minimum voltage loss in the contacts 
[4l] |. This assumption will allow us to switch to Dirich- 
let boundary conditions imposed on the potential profile, 
which ends up being computationally easier to handle. 
Since the contact potentials are assumed specified here- 
after, wc will use the MOM instead of explicitly solving 
Poisson's equation on a real space grid. 



IV. POTENTIAL SCREENING INSIDE A WIRE: 
THICK VS. THIN 

We will now concentrate on the potential profile in- 
side the wire, and discuss how its overall shape is deter- 
mined by broad, macroscopic geometrical features such 
as the wire thickness and screening length. We will see 
that screening is ineffective in a wire that is thinner than 
its Debye length. Later on, we will talk about finer de- 
tails, determined by the internal molecular structure of 
the wire itself. 
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FIG. 3: Electrostatic potential profile under two volts ap- 
plied bias for ballistic thin and thick wires (devices simulated 
shown as insets) using (a) semiclassical modeling with a con- 
tinuum structure; (b) quantum tight-binding modeling with 
an atomistic structure, and (c) fully ab-initio DFT modeling 
(adapted with permission from 0). Symbols represent the 
positions of the atoms along the device. The electrostatic po- 
tential drop along the wire decreases from thin to thick wire 
due to increased feasibility of increased transverse screening 
when the wire thickness W exceeds the Debye length \o, for 
which the field lines are shown schematically in (d). In the 
quantum calculations (b) and (c), the potential profile ex- 
hibits Friedel oscillations superposed on the long-wavelength 
Debye screening. Self-consistent charging under bias leads 
to a voltage asymmetry (see explanation in text as well as 
caption of Fig. 0J , although the I-V remains symmetric with 
respect to bias direction. 

Fig. ISIa) shows the potential profiles for a wire with 
a 2 volt applied bias, using a semiclassical, continuum 
approximation (Eq. The wire is 9 atoms long, while 
the contact block is 8 atoms long (only 3 atoms plotted) 
and 5x5 atoms in cross-section. The solid line shows the 
average electrochemical potential along the wire, which 
is spatially unvarying along the ballistic wire and drops 
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only at the contact interfaces. Using this electrochem- 
ical potential as a driving term in Eq. Ilbl gives the 
electrostatic potential profile </>. The Debye lengths of 
the contact and the wire regions are each assumed to be 
half the interatomic separation. For a thick slab repre- 
senting the wire (3x3 atoms in cross-section, geometry 
shown lower right), the potential profile (diamonds) is 
well-screened by charge rearrangement in the material of 
the wire described by the Poisson part (/)p (p) . As we pro- 
gressively thin the wire, however, we reach a point where 
the thickness of the wire is no longer much smaller than 
the Debye length (1x1 atom in cross-section, inset top 
left), at which point transverse charge screening become 
ineffective. This causes the field lines to penetrate in the 
transverse direction, generating a large ramp-like elec- 
trostatic potential profile along the wire (dashed line), 
essentially exhibiting just the Laplace solution (pL. A 
continuum analysis by Nitzan et al. leads to a similar 
conclusion, showing the dependence of the slope of the 
potential on the wire length, thickness and screening pa- 
rameter 0). 




E(eV) 

FIG. 4: Equilibrium transmission of thin and thick wires both 
show fine structure associated with Fabry-Perot end effects 
(blue solid line) superposed on the overall structure associ- 
ated with transverse modes (mode transmissions in red). The 
multiplicity of channels for a thick wire (subband structure 
shown in inset) arising due to quantization of the transverse 
(fcx) quasimomentum leads to the overall structure in its 
transmission (red). Under a small positive source bias, the 
drain electrochemical potential adds less negative charge to 
the thick wire than the source depletes, owing to the proxim- 
ity of the Fermi energy with a transmission mode bandedge 
that lies just above it in energy. The net positive charge floats 
down the molecular levels relative to Ef, giving rise to the 
asymmetry in the potential profiles in Fig. 

In Fig. |S{b), we arrive at the same conclusions, albeit 
from a more rigorous atomistic model Hamiltonian, in- 
cluding self-consisting charging effects and the detailed 



bonding geometry at the Au(lll) contacts, and employ- 
ing the full machinery of the quantum kinetic NEGF 
transport formalism fEas. [5I13I14|) . The same two cases 
are studied as before, with a thin wire constructed out 
of a chain of eight gold atoms connected to a triangle 
of FCC Au(lll) contact surface atoms (top left), and a 
thicker version generated by enclosing the central wire 
in a sheath of six identical gold wires around it (bot- 
tom right). The wires have quantum mechanical cou- 
plings in the transverse direction as well, allowing charge 
flow in that direction. Contacting the wires to the 3-D 
gold leads leads to a dimensionality mismatch and an as- 
sociated work-function mismatch, transferring about 0.8 
electrons at equilibrium to the thin wire and 2 electrons 
to the thick wire (the charging energy per electron L^o 
is about 2-3 eV for the thin and 1.1 eV for the thick 
wire, so the work function mismatch leads to a band ad- 
justment of about 2 eV for both thin and thick wires at 
equilibrium). The electrostatic potential profile along the 
central wire progressively changes from flat (diamonds) 
to ramp (dashed line) as the surrounding wires are re- 
moved |4lj[. A DFT version with a six atom gold wire 
(Fig. |3|;) |a| shows a similar profile. The screening along 
the transverse direction (Fig. |3|1) is feasible only if the 
wire thickness W exceeds the Debye length Ac- The os- 
cillations in the potential profile arise due to coherent 
Friedel oscillations (discussed later). A part of the po- 
tential drops across the end triangle of contact atoms 
(symbols 1 and 8 referring to the end wire atoms in b). 
Interestingly, even the surface potential along the thick 
wire (circles) shows screening although it is exposed on 
least one side. We believe the other wires screen the field 
over a large enough angle of the cross section that the 
field penetrating from the exposed part is minimal. 

The potential profile for the thin wire is spatially sym- 
metric, but that of the thick wire shows a significant 
amount of asymmetry. This asymmetry arises from self- 
consistent charging effects in the wire associated with its 
quantum capacitance. While the transmission of the thin 
wire is effectively structureless and approximately unity 
over a band (except for Fabry-Perot type resonances near 
the band-edges), the transmission of the thicker wire has 
considerable structure owing to the availability of var- 
ious transverse modes that kick in and out at various 
energies. Fig0| shows the zero bias, non self-consistent 
transmissions for a thin and a thick wire. The transmis- 
sions show a rapid oscillatory structure (blue solid lines) 
due to Fabry-Perot effects from the end superposed on a 
broader structure that arises from the transverse modes 
of the wires (red solid lines). For the thin 1-D wire, the 
transmission is unity between the band edges. For the 
thick wire, however, one can support multiple transverse 
modes associated with the seven wires (subband struc- 
ture shown in inset; there are seven subbands, of which 
some are doubly degenerate, denoted by 'x 2'). At 2 V 
bias the existence of an edge in the transmission mode 
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FIG. 5: Decrease in screening as a thick wire is pulled apart 
(same geometry and applied bias as in Fig. with only 
the part between atoms 1 and 8 shown here). For a closely 
packed coaxial wire, the potential profile of the central wire 
was screened (Fig. I^Jj, diamonds) compared to the thin wire 
(blue dotted line above, and blue dashed line in Fig. 03) • The 
screening is reduced drastically on pulling the wires just a lit- 
tle bit (left), since it eliminates the quantum couplings among 
the wires and prohibits rearrangement of screening charges in 
the transverse direction. The electrostatic coupling between 
wires is longer ranged, and is eliminated once the inter- wire 
separation D becomes comparable to the maximum distance 
between a charge and the nearest source/drain ground plane, 
i.e., half the wire length. The potential profile (right) of each 
wire now resembles that of an individual thin wire. 



spectrum above Ep causes the source electrochemical po- 
tential to add less charge to the wire than the drain re- 
moves, so about 0.08 electrons are removed from the thick 
wire. The average single-electron charging energy of the 
thick wire, ~ 1.1 eV/electron causes a shift in average 
potential profile associated with the charge depletion by 
~ 0.09 eV. The shift gets localized on only the central 
5 atoms or so where charge depletion is most, so the lo- 
cal electrostatic potentials at those points get lowered by 
about 0.14 — 0.2 V relative to Ep. as in Fig. mak- 
ing the overall electrostatic potential profile asymmetric. 
Such an asymmetry due to charging arises whenever the 
transmission spectrum has significant structure, as for a 
multimoded wire or a molecular conductor like PDT with 
distinct asymmetric highest occupied (HOMO) and low- 
est unoccupied (LUMO) levels. However, this charging 
effect is not readily observed experimentally, because it 
docs not affect the I-V directly. Although charging makes 
the potential profile asymmetric, the I-V is symmetric 
with drain bias because the sense of the asymmetry re- 
verses perfectly on reversing bias, so that the I-V does 
not show any sign of rectification ■ 

It is interesting to see how the screening effects get di- 
luted as the thick wire is progressively taken apart into its 
constituents. Fig. [3 (left) shows the hexagonal sheath of 
wires (cross-section shown in the inset) undergoing such 
a decimation process. First we lose the quantum cou- 
plings between adjacent wires, valid at a small increase 
in separation T) between them. This process in essence 
localizes the charges to the individual wires, disallowing 
any transverse charge rearrangement and making trans- 
verse screening quite inefficient. The corresponding elec- 
trostatic potential profile starts to look unscreened, like 
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FIG. 6: Dependence of the electrostatic potential on the di- 
mensionality of the electric field: (a) single atomic wire sim- 
ulated with 3-D Poisson's equation; (b) 2-D distribution of 
charges and fields, with no transverse variations, correspond- 
ing to a single SAM layer. Each wire atomic charge is smeared 
out into a line along the layer direction; (c) variation allowed 
only in 1-D, corresponding to a close-packed SAM or a molec- 
ular solid, with the charge on each atom smeared out into a 
2-D sheet along the SAM cross-section. A SAM screens fields 
in the transverse direction better, hence the trend in the po- 
tential profile with reducing dimensionality of the field lines. 



a ramp, being affected only by limited transverse flow 
of charge within each wire, around the individual atomic 
cores. The coupling at this stage is electrostatic, given by 
the inter-wire capacitances, which decrease linearly with 
T). Finally, when V k, L/2, where L is the wire length, 
the capacitive couplings with the source and drain elec- 
trodes dominate over the inter-wire capacitances, and the 
wire charges image on the nearer electrodes instead of on 
neighboring wires. At this stage, the assembly behaves 
in essence like unscreened wires, with each of their in- 
dividual potential profiles (Fig. O right) resembling the 
ramp-like profile of an isolated wire (blue dotted line). 

In this section, we saw that screening is ineffective for 
a wire that is thinner than its Debye length. This is be- 
cause of the feasibility of unscreened fields penetrating 
along the transverse direction (Fig. Old)) that can lead 
effectively to a poor overall screening. Transverse screen- 
ing can be improved by reducing the cross-section of the 
contacts. This eliminates field lines away from the con- 
tacts, generating a Laplace solution that is itself flat in 
between the two contacts. In our calculations, however, 
specifying the potential boundary conditions on the thin 
contact cross-sections alone is not sufficient to allow the 
method of moments to converge readily to a solution, 
due to incomplete specification of boundary conditions. 
Convergence requires fixing the potential at a large dis- 
tance, as if on imaginary gate electrodes. Similar conclu- 
sions have been obtained for ideal 1-D wires 01 • H has 
been shown, for instance, that long-ranged longitudinal 
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screening is ineffective in conductors such as metallic nan- 
otubes where the charge transport is in 1-D (along the 
axial direction of the nanotube surface). Thus doping 
the nanotube at an interface leads to an asymptotically 
nonvanishing surface charge along the nanotube |l7l |. 



pears well-screened |5j- 

VI. INTRAMOLECULAR POTENTIAL 
VARIATIONS 



V. POTENTIAL SCREENING BY THE 
ENVIRONMENT: SINGLE MOLECULE VS. 
MONOLAYER 

We have seen that the efficiency of voltage screening 
depends on the availability of adequate wire thickness to 
enable sufficient charge reorganization in the transverse 
direction. Thus for thin wires in isolation, such as in a 
break-junction experiment, screening is incomplete and 
the potential profile is essentially given by the Laplace 
part (f>L- In contrast, embedding the wire inside a self- 
assembled monolayer (SAM) allows sufRcient charge re- 
arrangement among wires, so that the transverse fields 
can be efhciently screened out by neighboring wires. We 
investigated transverse screening by neighboring wires by 
modeling a SAM that is very densely packed. For wire 
separation much smaller than the wire length (half the 
length to be precise, as we saw in the previous section), 
one could ignore charge and potential variations in the 
transverse direction. This process amounts to replac- 
ing the 3-D Poisson equation with a 2-D version that 
smears each gold atom into a transverse wire generat- 
ing a mat of gold atoms, or a 1-D version smearing each 
gold atom into a sheet generating a block or SAM of 
gold atoms (Fig. EJ. Instead of actually solving Pois- 
son's equation in reduced dimensions for each case, we 
replace the PPP 3-D Green's function in Eq.^jwith a 2- 
D version U{fi, rj) oc In |r^; — rj| and with the 1-D version 
U {xi, Xj) oc min(\xi — x j\, L ~ \xi — x j\) . For transport, 
we assume that the wires are separated and do not have 
strong transverse couplings, so that one does not obtain 
a transverse band of energies. Charge cannot flow easily 
in the transverse directions by hopping from wire to wire, 
and one gets essentially several copies of the same wire 
along the transverse directions. This means that instead 
of using _FiD and .Pjd functions along with the 2-D and 
1-D Poisson equations like we did in section IV, we use 
f{E)Ni and f{E)NA, where Ni and Na respectively de- 
note the number of wires per unit length and the number 
per unit area. 

The resulting self-consistent potentials arc shown in 
Fig. ini The inter-wire separations for the 2-D and 3-D 
cases are assumed to be equal to the separation between 
gold atoms (2.885 A) within the wire. From the figure, 
we see that as the Poisson equation is shifted to lower 
dimensions corresponding to the inclusion of more and 
more neighboring wires, the potential profile gets pro- 
gressively flatter due to enhanced screening efhcicncy of 
the transverse modes, so that the 1-D Poisson profile ap- 



The previous sections dealt with the broad shape of 
the potential profile across a molecular conductor. As 
was evident from Fig. |3 the qualitative features of the 
potential profile can basically be understood in terms 
of simple semiclassical continuum pictures, except for 
Fricdel oscillations and issues relating to the molecular 
density of states and single-electron charging. The ques- 
tion may arise at this stage as to whether the sophis- 
ticated machinery of quantum transport and NEGF, or 
the quantum Hamiltonian and atomicity of the device 
are important at all and whether they could have ob- 
servable consequences beyond what is predicted from a 
continuum description. The answer is that it could mat- 
ter, specifically for processes that involve the chemistry 
of the molecule. Specifically, we will discuss two such ex- 
amples: (i) potential barriers in PDT, and (ii) features 
related to the alignment of levels localized on different 
parts of the molecule, such as those generating a diode- 
like I-V or a negative differential resistance (NDR). 

Fig.dshows the electronic charge distribution on PDT 
at equilibrium and under a 2 volt bias. Fig. |H1 shows 
the corresponding electrostatic potential profile, includ- 
ing gold clusters at the end. The result, obtained from a 
semiempirical calculation, compares well with DFT cal- 
culations performed elsewhere 0. The largest voltage 
drop occurs between the end sulphur and gold atoms, 
with the sulphur atoms acting as a barrier (Fig. and 
Fig.lHJ, and minimal voltage drop within the gold cluster 
itself. The sulphur barrier is formed because of the ionic 
Au-S bonds, which involve transfer of electrons from the 
gold leads onto the electronegative sulphur atoms. Al- 
though it is hard to measure such an atomic barrier ex- 
perimentally, it is instructive to note that the barriers 
on sulphur do control the magnitude of the current to 
some extent, although the thinness of the barrier (few 
Angstroms) probably allows substantial electron tunnel- 
ing through it. Furthermore, the barrier regulates the 
flow of current, causing charge to pile up there under 
bias forming a residual resistivity dipolc, as the current 
rearranges to go around the benzene ring (Fig. (T)). 

Atomistic barriers and defects can regulate the direc- 
tionality and magnitude of current flow patterns within 
the molecule. In addition, coherent transport of elec- 
trons leads to Friedel oscillations in the potential pro- 
file around the defects and interfaces (FiglJb; also see 
|44|). This is a purely quantum mechanical phe- 
nomenon arising from the sharpness of the Fermi sur- 
face at low temperature. The Friedel oscillations are 
superposed on the long-wavelength Debye screening de- 
scribed earlier. The coherent oscillations can be elim- 
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(a) (b) 

FIG. 7: (a) Electronic charge distribution in PDT at equilibrium, smoothened artificially using a single Gaussian function for 
clarity. There is a distinct barrier on the sulphur atoms due to electron transfer from gold electrodes, (b) Change in electron 
distribution under bias shows an overall electron flow to the right. Charge piles up at points A and B generating two residual 
resistivity dipoles, as current needs to flow around the benzene ring. 




FIG. 8: (Calculated electrostatic potential profile across PDT 
connected to Au(lll) clusters shows maximum voltage drop 
across the sulphur barriers^with a linear Laplace potential (a 
DFT version is shown in [g). 

inated by incorporating incoherent scattering into the 
molecule through Biittiker probes 

The intramolecular electrostatic potential determines 
the local energy levels and the local density of states 
(LDOS), so that any substantial potential variation could 
bring various parts of the molecule in and out of res- 
onance. Such resonances can lead to observable conse- 
quences, such as generating an asymmetric, diode-like I- 
V by bringing the donor and acceptor levels at two ends of 
a molecule in and out of resonance (Aviram-Ratner mech- 
anism) analogous process is shown in Fig. |51 
where an artificially weakened bond in a QPC leads to 
maximum potential drop across it, separating the wire 
into two parts that separately equilibrate with the gold 
leads they are in contact with. The LDOS of the two 
parts slide past each other due to (f)L{r), bringing trans- 
mission peaks on both sides in and out of resonance. Such 




Weak Bond 

FIG. 9: Most of the applied bias drops across a weakened 
bond in a QPC (top left), effectively separating the two sides 
which individually equilibrate with the corresponding metal 
contacts. The QPC I-V shows a weak NDR (top right) due 
to voltage- induced alignment and misalignment of the LDOS 
on both sides of the defect (a DFT version is shown in [Tl|'l. 

a resonance leads to a weak NDR l33| (a DFT-based ver- 
sion was demonstrated earlier, in [4^). 



VII. INFLUENCE OF REMOTE CONTACTS: 
GATE MODULATION 

We have seen that the conductance of a molecule can 
be influenced simply by the Laplace part of the elec- 
trostatic potential, which depends on the geometry and 
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FIG. 10: Influence of a gate on the molecular I-V for (a) a gold wire and (b) a PDT molecule placed near the gate (gate is 
placed 1.5 nm away from the gold wire and 2.9 A from PDT). The gate controls the device potential and shifts the molecular 
levels relative to the contact electrochemical potentials. The current is unaffected for a gold wire due to its relatively featureless 
density of states that make it insensitive to level shifting. For PDT, the current tends to saturate when the drain electrochemical 
potential ventures into the HOMO-LUMO gap. The gate control, however, is quite poor due to the poor aspect-ratio in the 
problem, as well as metal-induced gap states (MIGS) from the gold contacts. Furthermore, the system is expected to have 
considerable gate leakage due to the thinness of the gate insulator (vacuum in this treatment), making the molecule a poor 
transistor. The gating can be improved in principle by using longer molecules with doped silicon contacts as source and drain, 
and high-k gate insulators to eliminate gate leakage currents. 



shape of the source and drain contacts. In addition, the 
potential profile in a conductor can be substantially influ- 
enced by the presence of a third (gate) electrode (Fig. ^1 
top left). Image charges on the gate electrode deempha- 
size the Poisson solution, so that the gate-induced molec- 
ular potential is essentially a Laplace contribution (jy^. A 
negative gate bias raises the molecular energy levels by 
increasing the average electrostatic potential (j)L (Fig. ^1 
top center) relative to the contact Fermi energy. For a 
molecule with constant DOS such as a thin gold wire, 
this shift does not affect the I-V, which remains ohmic 
with a quantized conductance Go = 2e^//i « 77/iS for 
varying gate voltage values (Fig. ^| top right). The rel- 
ative insensitivity of the conductance quantization to the 
contact geometry has been discussed at length in sev- 
eral papers [45l |. However, for a transistor involving a 
molecule such as PDT (Fig. bottom left) having a 
DOS with a lot of structure associated with broadened 
HOMO and LUMO levels, the overall potential shift af- 
fects the zero-bias conductance. For PDT, the conduc- 
tance increases with increasing negative gate voltage be- 
cause the molecule is essentially p-type (closest conduct- 
ing level to Ep is HOMO-based). In addition, the poten- 
tial profile is skewed towards the source-end (Fig.^| bot- 
tom center), which has a common ground with the gate 
electrode. Good gate control in a "well-tempered" metal 
oxide semiconductor field effect transistor (MOSFET) 
makes the Laplace potential insensitive to source-drain 
bias and sets it by the source-gate bias instead (as in 
Fig.lTUIl [iif. The asymmetric potential profile effectively 



tends to pin the channel potential to the source electro- 
chemical potential, so that only the electrochemical po- 
tential of the drain electrode varies under source-drain 
bias. For negative drain bias, the drain electrochemical 
potential iijj enters the molecular HOMO-LUMO gap 
(HLG), leading to an I-V characteristic with decreasing 
slope due to a decreasing DOS in the gap (Fig.^1 bottom 
right). Different gate voltages yield different zero-bias 
values Ep for the drain electrochemical potential rel- 
ative to the levels. This gives a gate- voltage dependence 
of the saturation current, leading to MOSFET-like I-V 
characteristics with gate modulation ^3 ■ 

While the above gate control mechanism describes the 
principle of operation of an ideal ballistic silicon MOS- 
FET |4g, the quantitative conclusions (saturation, gate 
modulation) are usually severely compromised when the 
MOSFET is scaled to molecular dimensions. Electro- 
static gate control requires the gate electrode to be 
substantially closer to the molecule than the source- 
drain electrodes, while good saturation in the IVs (large 
impedances) require a vanishing DOS in the HLG. For 
long molecules [l3 one can have a modest oxide thick- 
nesses that still yields appreciable gate control. However, 
for small molecules ~ 10 A, the oxide thickness cannot 
be scaled down enough without causing dielectric break- 
down, leading to poor electrostatic gate control Fur- 
thermore, the current saturation is poor (Fie. 110(1 due to 
a non-negligible gap DOS arising from the broad HOMO 
tails generated by MIGS from the gold contacts. Improv- 
ing the gate control and the impedance require the uti- 
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lization of high-k gate dielectrics and degenerately doped 
semiconducting contacts. Alternately, one could envi- 
sion transistor action based on non-electrostatic princi- 
ples, such as by gating the molecular conformations, for 
example |50| . 

The gate control can be enhanced by wrapping the 
molecule with cylindrical gates surrounding the wire. 
Such a scheme is employed in silicon transistors to pro- 
duce dual-gate MOSFETS FINFETs ^ or struc- 
tures with wrap-around gates [s^- We have modeled 
the effect of multiple gates by placing two, three and 
four rectangular gate electrodes symmetrically around 
the wire. Increasing the number of gate electrodes in- 
creases the potential barrier, along with an overall in- 
crease in the average potential. This leads to superior 
gate control. The practicality of this scheme, however, 
depends on the capacity to have multiple gates closely 
spaced around the molecular wire and insulated from it. 



IX. CONCLUSIONS 

The electrostatic potential profile across a conductor 
controls the distribution and magnitude of current across 
it. While broader features such as the overall conductiv- 
ity and screening lengths determine the shape of the I-V 
curves in many cases, atomistic features involving level 
resonances can influence these I-Vs in a nontrivial way, 
necessitating a proper atomistic, quantum kinetic treat- 
ment of charge transport in molecular devices. We have 
performed such a treatment within a coupled Poisson- 
NEGF formalism, while simplifying the molecular chem- 
istry to essentially one orbital per atom for computa- 
tional simplicity and ease of illustration (a more sophis- 
ticated DFT version for some of these profiles has been 
demonstrated elsewhere) . We have seen that the Poisson 
part of the electrostatic potential profile depends on the 
strength of transverse and longitudinal screening, and 
can yield nontrivial charging effects. A gate electrode 
generates image charges that deemphasize the role of 
the Poisson part, while modulating the Laplace solution 
suitably to significantly alter the device TV characteris- 
tic and yield saturating currents. However, feasibility of 
electrostatic gate control depends on the channel length 
to oxide thickness aspect-ratio in the problem, and may 
need to be supplanted by alternate physical principles, 
such as employing conformational degrees of freedom for 
example. 
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